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ABSTRACT 

Analysis of a longitudinal wave event observed by the Atmospheric Imaging Assembly (AIA) onboard 
the Solar Dynamics Observatory (SDO) is presented. A time sequence of 131 A images reveals 
that a C-class flare occurred at one footpoint of a large loop and triggered an intensity disturbance 
(enhancement) propagating along it. The spatial features and temporal evolution suggest that a 
fundamental standing slow-mode wave could be set up quickly after meeting of two initial disturbances 
from the opposite footpoints. The oscillations have a period of ^12 min and a decay time of ^9 min. 

The measured phase speed of 500T50 km s“^ matches the sound speed in the heated loop of ^10 MK, 
confirming that the observed waves are of slow mode. We derive the time-dependent temperature and 
electron density wave signals from six AIA extreme-ultraviolet (EUV) channels, and find that they 
are nearly in phase. The measured polytropic index from the temperature and density perturbations 
is 1.64 ± 0.08 close to the adiabatic index of 5/3 for an ideal monatomic gas. The interpretation 
based on a ID linear MHD model suggests that the thermal conductivity is suppressed by at least a 
factor of 3 in the hot flare loop at 9 MK and above. The viscosity coefficient is determined by coronal 
seismology from the observed wave when only considering the compressive viscosity dissipation. We 
find that to interpret the rapid wave damping, the classical compressive viscosity coefficient needs to 
be enhanced by a factor of 15 as the upper limit. 

Subject headings: Sun: Flares — Sun: corona — Sun: oscillations — waves — Sun: UV radiation 


1. INTRODUCTION 

Coronal seismology is a technique for measur¬ 
ing physical quantities of the corona by matching 
observations with magnetohydrodynam ic (MHD) the¬ 
ory of waves in structured plasma (jRpberts et al. I 
Il983f ). Considerable progress has been made in 
the study of coronal seismology in the last d ecade 
(see reviews bv [f^kariakov fc Verwichte I I2QQ5I : 
iDe Moortel fc Nakariak^ ] l 2 nn ILiii fc OfmaiTI ImTl . 
From o bservations of propagatin g slow magnetoacoustic 
waves, IVan Doorsselaere et al. H (l2nTl estimated the 
polytropic index to be ^1 in the warm corona. The 
knowledge of the appropriate value of the polytropic 
index is important in hydrodynamic and MHD models 
of the solar and stellar coronae a s well as of space 
plasmas (e.g. iPnc ovkin et al. I Il997l : iRilev et al. I I2QQ1I: 
iJacobs fc Poedts~ll2Qll[). 

Longitudinal hot loop oscillations were first discov¬ 
ered with the SUMER/S'OiLO as periodic variations of 
Doppler shift in F exix and Fexxi lines (jWang et al. I 
I2QQ2I : IWang I |2Q11[ ) . These oscillations were mainly in¬ 
terpreted as fundamental standing slow modes because 
their phase speed is close to sound speed in the loop, 
and there is a quarter-period phase shift between the ve- 
l ocity and intensity oscillations in some observed cases 
(jWang et al. ll2QQ3allbl ). The initiation of the waves was 
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often associated wi th small flares at the loop footpoint 
(jWang et al. ll2QQ^ . 3D MHD simulations show that the 
fundamental standing slow mode wave can be excited by 
a velocity pulse or impulsive onset of flows near one foot¬ 
point of the loop (jSelwa fc Ofman I I2QQ9I: lOfman et al. I 
l2Q12f ). These hot loop oscill ations typically show a rapid 
decay. MHD simulations bv lOfman fc Wan^ (j2QQ2f ) sug¬ 
gested that thermal conduction is the dominant wave 
damping mechanism enhanced by nonlinear effect. The 
role of other physical effects such as compressive viscos¬ 
ity, non-equilibrium ionization, shock dissipation, loop 
co oling, etc was also studied theoretically (see a review 

bv lWang 


11^ IR 

Iv. iKnms 


Recently, iKnmar et al 


.nderman ll2Q13t lAl-Ghafri et al. II2Q141). 

120131 ) reported a longitudinal 


wave event observed with the SDO/AIA, showing simi¬ 
lar physical p roperties as found p reviously in SUMER 
observations (jWang et al. II2QQ3H ). However, the AIA 
wave was seen bouncing back and forth in the heated 
loop, suggesting that it may be a propagating mode in 
contrast to the standing modes identified by SUMER. In 
this study we report the first SDO/AIA case that shows 
clear signatures in agreement with a fundamental stand¬ 
ing slow mode wave. We find that the temperature and 
density perturbations in >9 MK plasma are nearly in 


phase and the measured polytropic index accords well 
with the classical value (5/3) of the adiabatic index for an 
ideal monoatomic gas, suggesting that the thermal con¬ 
ductivity in hot plasma is much weaker than predicted 
by the classical theory. This new finding may challenge 
our current understandin g of thermal energy transpo rt 
in solar and stellar flares (jShibata fc Yokovama l[2QQ^ . 


2. OBSERVATIONS AND DATA ANALYSIS 
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Fig. 1.— Observations of a longitudinal wave event on 2013 December 28. (a) AIA 131 A image. A sampled slice (green lines) along the 
oscillating loop is used to construct the time-distance plot shown in Figure [2] The white thick curve indicates the reconstructed 3D loop 
(same as the sampled track) and the dashed curve for a fitted circular model, (b) AIA 131 A base-difference image relative to a pre-flare 
image at 12:40 UT. The inset shows a co-temporal Hinode/XRT image in the Be_med filter, (c) and (d): The detrended 131 A images at 
12:51 UT and 12:57 UT, showing the antiphase intensity perturbations at the opposite legs of the loop. (Animations are available in the 
online journal) 


A GOES C3.0-class flare occurred on 2013 Decem¬ 
ber 28 in NOAA Active Region (AR) 11936 near disk 
center of the Sun. The flare began at 12:40 UT and 
peaked at 12:47 UT. This flare excited longitudinal waves 
in a heated loop observed with 5'DO/AIA. The AIA 
records the full Sun in ten extreme ultraviolet (EUV) and 
UV passbands with high spatio-temporal resolution and 
wide temperature coverage flogT=3.7-7.3: iLemen et al. I 
l2012[) . We analyze the loop dynamics and thermal prop¬ 
erty using the AIA data. 

The longitudinal waves were clearly detected in AIA 
131 A (dominated by Eexxi, formed at ^11 MK) band. 
The animation shows that a large loop that connected the 
flare site at one end (marked with a curved cut in Eig- 
urema)) developed quickly and exhibited intensity oscil¬ 
lations along the loop. Similar emission features seen in 
the AIA 131 A base-difference image and a co-temporal 
soft X-ray image (F igure HTb)) observed by Hinode/XRT 
(|Golnb et al. ll2QQ^ indicate that the oscillating loop is 
hot. To emphasize the evolution of intensity perturba¬ 
tions, we processed the 131 A images by subtracting the 


slowly-varying trend at each pixels, and then normal¬ 
ized the perturbation to the trend, where the trend was 
derived using the Fourier low-pass (>20 min) filtering. 
The detrended images show alternate intensity enhance¬ 
ments at the opposite legs with rapid decay (Figures [TJc) 
and (d)). The animation shows that the loop oscillations 
lasted for about two periods before they faded out. 

To measure the wave properties we sampled the 131 A 
images from a 21-pixel (13 arcsec) wide slice (outlined in 
Figure [TJa)) over a period of 12:40-13:20 UT. We aver¬ 
aged the slices over the width and stacked them in time to 
construct a time-distance plot (Figure [2j a)). Figure [2]^b) 
shows the corresponding detrended time-distance plot. 
The intensity oscillations along each leg are nearly in- 
phase but anti-phase between the opposite legs match¬ 
ing the signature for a fundament al standing slow mode 
predi cted by forward modelings (|Tarovan fc Bradshaw I 
|200l). Figure m; c) shows the time profile of mean in¬ 
tensity for a cut selected at the leg with the brightest 
emission (at segment 15 shown in Figure [3fa)). We fit- 
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Fig. 2. — Measurements of the loop oscillations seen in 131 A. (a) Time-distance plot for the sample wave path traced out in Figure fTlah 
showing wave perturbations along the loop. In the plot, the distance is measured from the end near the flare site, (b) Same as (a) but with 
the slowly-varying trend subtracted at each position. The two arrows marked indicate two initial brightenings propagating in the opposite 
directions, (c) Time profile of the average counts extracted from a 9-arcsec wide region (narrow box in (a)). The solid curve shows a best 
fit which includes a parabolic trend (dashed line), (d) Detrended light curve with the best fit (solid line). 


ted the oscillatory signals with the function 

I{t) = Imsin exp + Io{t), 

( 1 ) 

where Im^ i", 0, and to are the amplitude, period, 
decay time, initial phase, and reference time, respec¬ 
tively, and /o(t) is a parabolic trend. The best fits to 
the light curve with and without the trend are shown in 
Figures [2fc) and[2fd), respectively. The period of oscil¬ 
lation is P = 12.0 ±0.1 min, the decay time r = 9.7±0.3 
min, and the initial amplitude relative to the trend is 
/m//o(to) =0.65 ±0.01. 

The phase speed is important for identification of 
wave modes, which requires an estimate of the loop 
length. We estimated the loop 3D geometry using the 
curvature radius maximization method which assumes 
the line-of-sight (LOS) coordinates of the observed loop 
to be the corresponding LOS coordinates of a circn - 
lar model (|Aschwanden et al71l2QQ2l: lAschwanden 1120091) . 
Figure [TJa) shows the fitted circular model and the so¬ 
lution of the 3D geometry which has an identical 2D 
projection as observed. We obtained the loop length 
I/=179 Mm, about 33 Mm longer than that of its 2D 
projection. The wave phase speed was estimated by 
Vp=2LlP = 500 ± 50 km s“^, where an uncertainty in L 
is taken as O.IP. The phase speed is close to the sound 
speed of 480 km s“^ for the loop temperature T ^10 MK, 

estimated using Cs = 152Tg^^ km s“^, where Tq = T/[l 
MK] and taking the adiabatic index 7=5/3. This result 
supports the interpretation of the observed waves as a 
fundamental slow mode. 

To conduct accurate diagnostics of the electron tem¬ 
perature and density of the oscillating loop we utilized a 
regularized differential emission measure (DEM) analysis 
on AIA images in six EUV channels (94, 131, 171, 193, 
211, and 335 A) ([Hannah fc Kont^ 120121. I2013D . As 
an illustration, we chose the images at ^12:53 UT when 
the entire loop became obvious to trace, and divided the 
entire loop into 23 subsections (Figure [3f a)). We per¬ 
formed the DEM inversion for each segment and fitted 
the DEM curve with a triple-Gaussian model to sepa¬ 


rate a hot component {^hot{T)) from the background by 
assuming that the hot component comes from the fore¬ 
ground hot oscillating loop (Figures [3l^b) an d (h)). The 
applic ation of this technique was detailed in iSnn et alTI 
(|2Q13f ). Figures ^d) and (e) show the loop tempera¬ 
ture (the centroid of the hottest Gaussian component) 
and the loop emission measure {EM = f ^hot(7")dT, the 
pink area in Figure [3l^b)). We measured the loop width 
by fitting the cross-sectional flux profile with a Gaus¬ 
sian function, and obtained the mean FWHM width, 
w = 13.8 ± 1.7 Mm for segments 8-18 (Figure [3l^f)). By 
using n = yjEM/w with a filling factor of unity (lead¬ 
ing to a lower limit of n), the loop electron density was 
estimated (Figure [3l^g)). Both the loop temperature and 
density are found to vary in a small range along the loop 
(T=7-12 MK and n=(1.6-3.0)x 10^ cm-^). 

To study the time evolution of thermal properties, we 
focused on segment 15 with the maximum EM and used 
its average flux in six AIA bands to perform a time series 
of DEM inversions with a cadence of 24 s. By applying 
the triple-Gaussian DEM analysis (Figures [3fb) and (c) 
and the animation in the online version), we obtained 
time profiles of the loop temperature and electron den¬ 
sity (Figures m^a) and (b)), where the loop width w was 
taken as a constant (^14 Mm) in calculation of the den¬ 
sity evolution. The mean temperature over the lifetime 
of oscillations is found tobeT = 8.7±1.7 MK, and the 
mean density n = (2.6 ± 0.2) x 10^ cm“^. The wave sig¬ 
nals in measured temperature and density are evident. 
We fitted the oscillations with a damped sine-function 
of the same form as equation O- The fitted oscillations 
and the background trends are shown in Figures SJ^a) and 
(b). The detrended oscillations with the fitted parame¬ 
ters are shown in Figures EJ^c) and (d). We find that the 
temperature and density oscillations have similar periods 
and they are nearly in phase. The initial phase shift is 
A0 = 6° ±23°. The phase shift measured using the cross 
correlation from the relative perturbations to the trend 
is about 12°. 

Under the poly tropic assumption p ^ with a poly¬ 

tropic index a, it can be derived from linearized ideal 
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Fig. 3.— The DEM analysis of loop thermal properties, (a) AIA 131 A subimage of the analyzed loop. The traced loop is divided into 
23 segments (solid lines) with segment z=0 near the flare site. Each segment has a size of 15x21 pixels. A 70-pixels wide region outlined 
with dashed lines is selected for measuring the width of loop cross sections. The loop segment z=15 is used for the analysis of EM and 
temperature evolution shown in Fig. [4] Panels (b) and (c) show the regularized DEM inversion result for segment 15 at 12:53 UT and 17 
min later. The DEM profile (crosses) is fitted to a triple-Gaussian model (curves) for isolating the foreground hot loop contribution (pink), 
(d) Loop temperature. The gray band shows the temperature range (Tcrx-, Gaussian width in LogT). (e) Loop EM. (f) Loop width, (g) 
Electron density, (h) Reduced chi-squares (x^) for the regularized DEM and the Gaussian model in data space. Error bars in (d)-(g) are 
the 1(7 fitting error. (An animation is available in the online journal) 

MHD theory that (e.g. I Van Doorsselaere et al. ll2Qll[) : based on the wave amplitude measurements as 


t' 


{a 



( 2 ) 


where p, n, T are the gas pressure, electron number den¬ 
sity, and temperature, respectively. A superscript ( ) in¬ 
dicates perturbed quantities and a subscript (o) stands 
for equilibrium quantities. The process is adiabatic only 
if ^=7=5/3. The polytropic index can be measured us¬ 
ing the linear relationship between the observables T /Tq 
and n I no. Note that the equation (j2j) does not generally 
hold true in the nonideal MHD case. For example, ther¬ 
mal conduction can lead a large phase shift between T 
and n (see the discussion in Section [3j). In this case, the 
polytropic index should be measured either using equa¬ 
tion © for the data after removing the phase shift, or 


^ AT{t) 

An{t) 


+ I5 


(3) 


where and An{t) are the decaying temperature 

and density amplitudes normalized to the correspond¬ 
ing trend. Here we chose equation m to measure the 
polytropic index because the observed temperature and 
density perturbations are nearly in phase and the ad¬ 
vantage compared to using equation is that the as¬ 
sumption of the oscillations as a damped sine-function 
is not required. We obtained a = 1.64 ± 0.08 using the 
linear least-squares fitting to the scatter plot of the data 
(T /To{t) and N /No{t)) over the wave lifetime of 12:46- 
13:20 UT (Figure [5ja)), and obtained a = 1.66 ±0.09 by 
fitting only the data within the first oscillation period 
(Figure mb)). We find that the value of a agrees well 
with the adiabatic index of 5/3 for fully ionized coronal 
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plasmas. Note that the assumptions of the line-of-sight 
column depth and the filling factor used in determina¬ 
tion of the loop density from EM have little effect on the 
accuracy of the measurement because a depends only on 
the relative perturbations of the density (n /no). 

3. DISCUSSION AND CONCLUSIONS 

We have studied a non-eruptive ffare event (not as¬ 
sociated with a CME) occurring near one footpoint of 
a large coronal loop observed by SDO/AIA. The ffare 
impulsively heats the loop to above 10 MK and also gen¬ 
erates the longitudinal waves in the hot loop. Despite 
of difficulties in distinguishing propagating and stand¬ 
ing waves based on only intensity information, we ar¬ 
gue that the observed wave is most likely of a fun¬ 
damental standing slow mode. This is suggested by 
the signature that the temperature and density oscilla¬ 
tions match well a harmonic wave with the frequency 
close to that for the fundamental mode. In contrast, 
a reffecting single wave pulse does not follow a sine- 
function in temporal variati on and its time distance plot 
exhibits a ‘zigzag’ pattern (iDe Moortel fc Hood I I2QQ3I: 
iTarovan fc Bradshaw I|2QQ1 : iKumar et al. ll2Q13[ ). In ad¬ 
dition, the nearly inphase intensity oscillations shown 
along each leg after ^12:50 UT suggest that the stand¬ 
ing wave appears to form within one wave period after 
the ffare, consistent with the case of standing slow waves 
observed by SOHO/SUMER 
The quick setup of the standing wave could be related to 
meeting of the two oppositely-propagating brightenings 
marked with arrows in Eigure[2jb)) in the initial phase 
Selwa fc Ofman |[2QQ9[ ). however, a detailed study on the 
wave excitation is beyond the scope of this Letter. 

Here, we focus our analyses on the plasma thermal and 
wave properties of the hot loop, and have measured the 
polytropic index a from the temperature and electron 
density perturbations. We find that the value of a is close 
to 5/3 within a small uncertainty particularlly in the 
first oscillation period. This is striking because it implies 
that the energy equation in the MHD equations can be 
well represented with an adiabatic form, or the non-ideal 
effects such as energy gain and losses on the timescale 
of oscillations are negligible in the hot plasma of T >9 
MK. The reasons are detailed in the following based on 
the linearized ID hydrodynamic equations including the 
effects of gravity, heat conduction, compressive viscos- 
ity, and optically thin radiative losses and gain s (e.g., 
iDe Moortel fc Hood ll 2 QQ^ iSigalotti et al. 1I2QQ7I ). Note 
that both the measurements of a and the theoretical in¬ 
terpretations below are independent on whether the ob¬ 
served slow wave is of propagating or standing mode be¬ 
cause they are essentially identical in physical property. 

(1) We know that for an isothermal loop in hydrostatic 
equilibrium, the pressure scale height H = SOTe Mm. Eor 
the hot loop at T MK, we have H/h «9 where h=50 
Mm, the estimated loop height. Thus the stratification 
effect here can be neglected. 

( 2 ) We estimate the effect of radiative loss from the 
ratio of the oscillation period to the radiation time scale 

r = Here Trad = 345 Qn^^T^^^ s, where ng = 

n/[I0^cm“^] (|Snn et al. II2QI3[) . Eor the measured pa¬ 
rameters T=8.7 MK, n = 2.6 x 10^ cm“^ and P=I2 min, 
we obtain Trad ~570 min and r «0.02. Thus radiative 



cooling is negligible on the oscillation period timescale. 

(3) It is known that when thermal conduction domi¬ 
nates in the energy equation, it introduces a phase shift 
( A(^) between the density and temperature perturb ations 
(|Owen et al. ll2QQ^ IVan Doorsselaere et al. ir 2 QIl[) . The 
linear approximation gives the following relations 


tanA0 

At 


7r(7 - 1)^11 
kBC^Pn 

(7 — I)cosA0 


[= a - 1 ], 


(4) 

(5) 


where 7 = 5 / 3 , ks the Boltzmann constant, = 7.8 x 
10 - 72 ^ 5/2 cm“^s“^K“^ the thermal c onductivity 

parallel to the magnetic field (|SDitzer Ill962[) . Note that 
the above equations hold true no matter whether ther¬ 
mal conduction or compressive viscosity dominates in 
the wave damping (see item (4)). ID MHD simulations 
have shown the presence of a large phase shift between 
T and n in hot loop s due to thermal cond uction (e.g. 
lOfman fc Wang I [20021 : iSkalotti et al. I [ 200 ^ . 

However, we find that the temperature and density 
oscillations are nearly in phase. This result suggests 
that the thermal conduction may be much weaker than 
the classical theory predicted in the hot fiaring loop. 
With the loop parameters used in item (2) we obtain 
the theory-predicted A(/) = 35° using equation (j4)). The 
comparison of the observed phase shift (^ 12 °) with the 
predicted suggests a reduction of (parallel) thermal con¬ 
ductivity by a factor of three. Moreover, from the mea¬ 
sured polytropic index of a=1.64 we derive the phase 
shift to be ^ 16° using equation m, which is consis¬ 
tent with the directly measured value. Thus both the 
phase shift and polytropic index measurements support 
the conduction suppression in this event. In addition, the 
results of a=1.66 that is even closer to 7 = 5/3 and a bet¬ 
ter in-phase relationship between T and n within the 
first oscillation period suggest that thermal conduction 
tends to be more suppressed when the plasma is hotter 
(T >9 MK). 

In addition, we notice that for hot loop oscillations ob¬ 
served by SUMER and SDO/AI A (jWang et al. ll2QQ3i] . 
l2QQ7l : iKumar et al. I l2QI3l . l2Q15f ). the measured phase 
speeds are close to the adiabatic sound speed within un¬ 
certainties, also suggesting a ^5/3. This is in contrast to 
th e theoretical prediction from linear MHD simulations 
bv iDe Moortel fc Hood I (|2QQ3f ). They showed that when 
thermal conduction is dominant due to the high tempera¬ 
tures, the perturbations propagate largely undamped, at 
the slower, isothemal sound speed, which implies a ^1. 

(4) In using equation (j2j) to measure the polytropic 
index, we have assumed that the loop temperature (or 
sound speed) varies gradually with time relative to the 
wave period (known as the WKB approximation). To 
check its validity, we quantify the “gradual variation” 
as the requirement that changes in wave period must 
be small over a single period. The condition can be 
expressed as \dP/dt\ I. The fitting to the evolu¬ 
tion of the loop temperature (Eigure SJ^a)) gives the 
trend, T{t) = Tg + at + 6 t^, where the initial temper¬ 
ature To=12 .6 MK, the polynomial coefficients a=—0.38 
MK min“^ and 5=0.0064 MK min“^. Considering 

P(t) = 2T/cs(T(t)), we obtain \dP/dt\ = PoTQ^^^|0.5a + 
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Fig. 4. — Evolution of thermal properties for the loop segment 15. (a) Temperature. The red solid curve shows the best fit to the 
oscillatory signals, and the dashed curve is the parabolic trend, (b) Same as (a) but for electron density. The gray band and error bars in 
(a) and (b) have same meanings as in Figure [3] (c) Detrended time profile of the temperature (crosses) and the best fit (red curve), (d) 
Same as (c) but for electron density. 




Fig. 5.— Measurements of the polytropic index a. (a) The scatter plot of electron density and temperature perturbations (pluses), with 
the best-fitting line (solid) and the line of 7=5/3 (dotted). The dashed lines indicate the Tier fitting error, (b) Same as (a) but for the 
data within 12 min after 12:46:20 UT. 


btlT(t)-^/^, where Pq = 2L/cs(To)=10.5 min. The cal¬ 
culation of its maximum value |dP/dt|niax=0-16 indicates 
that the WKB approximation is valid in this case. This 
justifies the removal of the trend from the oscillations in 
our analysis. 

(5) As the viscous heating belongs to the second-order 
term in the linearized energy equation, its effect on wave 
damping can be neglected. Thus the only dominant 
damping mechanism left is the momentum loss by the 
viscous forces. From the dispersion relation derived from 
the velocity wave equation whe n only compressive vis - 
cosity is present (see Eq.(33) in iSigalotti et al.~l I2QQ7I ). 
we can obtain the coefficient of compressive viscosity in 
terms of the observables by 

_ SjpoT _ SjkBnoToT 

“ 87r2(r/P)2 +2 “ 47r2(T/P)2 + 1’ ^ ^ 

where 7=5/3 and po = 2nokBTo. By taking P = 12.0 ± 


0.6 min and r = 9.2 ±2.1 min (the mean for temperature 
and density oscillations) with Tq = 8.7±1.7 MK and no = 
(2.6 ± 0.2) X 10^ cm“^, we obtain rjs = 356 ± 195 g cm“^ 
s“^. In comparison, we calculate the classical Braginskii 

5/2 

compressive viscosity coefficient using rjo = O.lTg^ g 
cm“^ s“^ (e.g., lOfman fc Wang |[2002[) . and obtain 770 ~ 
23 g cm“^ s“^, thus, r]s/Vo = 15. This implies that to 
interpret the wave damping timescale by the compressive 
viscosity alone the classical viscosity coefficient needs to 
be enhanced by a factor of 15, which may be regarded as 
an upper limit considering the additional effects such as 
weak nonlinearity, thermal conduction and stratification. 

In summary, we have found quantitative evidence of 
thermal conduction suppression in a hot flare loop by 
coronal seismology of the slow-mode waves. This re¬ 
sult suggests that the flare loop should cool much slower 
than expected from the classical Spitzer conductive cool- 
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ing. Our studied event is indeed of such long-duration 
events (LDEs) which are the flares with a slower-than- 
expec ted decay rate in soft X-ray and EUV radiation 
lEorbes fc Acton 119961: iTakahashi fc Watanabe I 
l2Q12 b To explain the LDEs and 
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flare loop-top sources, some previous studies have sug¬ 
gested the mechanism of continuous he ating (jWarren I 
l2QQ6l iLiu et al~1 12Q13I: iSnn et aI71l2Q13l) or conduction 
suppressi on (|McTiernan et alTI 119931: 1 Jiang et al. 

iLi et al. 1 l2fll2D . Our study confirms the effect of the 

latter mechanism in a more direct way. Laboratory ex¬ 
periments and numerical studies showed that the ac¬ 
tual conductivity is smaller (by at least a factor of 
two) than that given by Spitzer when I < 30A where 
I is the temperature gradient scale le ngth and A the 
mean free path of thermal electrons (IRcll et, a,lJ 
iLuciani et al. We estimate Z/A ~ 30 in this case, 

suggesting that the nonlocal conducti on may account for 
the observed conduction suppression (jMatte fc Virmo^ 


Il982l : [Rosner et al. 1119861). Bv studying t he evolution of 
flare loop-top sources, [jiang et al. I (j2QQ6f ) suggested that 
plasma waves or turbulence may play an important role 
in suppressing the conduction during the decay phase of 
flares. The mechanism is similar to that used for inter¬ 
preting the significant reduction of thermal conductiv¬ 
ity in galaxy-cluster cooling flows by a tangled magnetic 
field (jChandran fc Cowlev 11199'^ . Einally, we conclude 
that the result of conduction suppression may also shed 
light on the cor onal heating problem (see the review by 
lKlimchuOl2Q15[) , because weak thermal conductivity im¬ 
plies smaller conductive losses, and an extended lifetime 
of individual nanoflares, increasing the average coronal 
temperature for the same heating rate. 
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NNX12AB34G and the NASA Cooperative Agreement 
NNGIIPLIOA to CUA. L. O. and E. P. acknowledge to 
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mission for NASAs Living With a Star (LWS) program. 
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